Thermalization and dissipation in out of equilibrium quantum systems: A 
perturbative renormalization group approach 
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A perturbative renormalization group approach is employed to study the effect of a periodic po- 
tential on a system of one-dimensional bosons in a non-equilibrium steady-state due to an initial 
interaction quench. The renormalization group flows are modified significantly from the well known 
equilibrium Berezinski-Kosterlitz-Thouless form. They show several new features such as, a gen- 
eration of an effective temperature, generation of dissipation, as well as a change in the location 
of the quantum critical point separating the weak coupling and strong coupling phases. Detailed 
results on the weak-coupling side of the phase diagram are presented, such as the renormalization 
of the parameters and the asymptotic behavior of the correlation functions. The physical origin of 
the generated temperature and friction is discussed. 
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I. INTRODUCTION 

The high degree of tunability and control associated 
with cold-atomic gasesi has motivated an explosion of 
theoretical activity involving the study of dynamics of 
interacting quantum systems. Several interesting prob- 
lems could be studied in these systems such as quantum 
quenches (see Ref. 12 and references therein) and other 
classes of steady-state nonequilibrium phenomena such as 
systems subjected to a time dependent noised In all these 
situations one of the fundamental questions is what is the 
steady-state of the system, and in particular whether the 
system can thermalize. 

Not surprisingly in such a difficult problem, exactly 
solvable models and integrable systems have proven to 
be a good playground to address these issues. The study 
of quench dynamics in many integrable systems reveal 
that the large number of conserved quantities in the sys- 
tem prevents the system from thermalizing. Instead the 
long time behavior of the system is characterized by time- 
dependent or time-independent nonequilibrium states4^ 
Often after a time-averaging, the resultant steady-state 
can be described by a generalized Gibbs ensemble (GGE) 
constructed from identifying the conserved quantities of 
the system^ - — However the generality and applicability 
of the GGE remains under debate since not all observ- 
ables can be described by itJ^r— 

Besides a lack of complete understanding of the long 
time behavior of integrable and even some exactly solv- 
able models after a quench, the more complex question 
of how the presence of non-trivial interactions that cause 
scattering and/or break integrability affect the behavior 
is largely open. Numerical studies on finite systems show 
that for large enough breaking of integrability, thermal- 
ization is associated with a change of the level statis- 
tics from Poisson to Wigner- Dyson,— with the onset of 
thermalization occurring via a two-step process, where 
in the first step the system is trapped for a long time in 



a nonequilibrium prethermalized state^2r— At the same 
time, fluctuations in a finite-size system also play an im- 
portant rol o 25 ' 26 so that generalization of numerical re- 
sults to systems in the thermodynamic limit may not be 
straightforward. 

The aim of the current paper is to address some of 
these questions in precisely this limit of long times and 
infinite system size where numerical studies are hard. In 
particular we study an exactly solvable model which is in 
a nonequilibrium steady-state due to an initial quench, 
and explore the stability of the resulting athermal state 
to non-trivial interactions that generate mode-coupling. 
The effects of mode-coupling will be treated within a per- 
turbative renormalization group (RG) approach. While 
there are many candidate models to study this physics, 
due to experimental relevance, and its relative simplic- 
ity, we choose to study a one dimensional system of in- 
teracting bosons, leading to the so called Luttinger liq- 
uid physics.2£ The excitations of such a system can be 
represented by density modes, which are essentially in- 
dependent. On such a system, quenches corresponding 
to a change of the interaction reveal a steady-state which 
still has independent modes, but these are now character- 
ized by a nonequilibrium distribution that does not relax 
to a thermal state ,H5i2§r.2i Note that we use the term 
"steady-state" to reflect the fact that averages of various 
physical observables reach a time-independent value at 
long times after the quench. 

We study the effect of a mode coupling term, such as 
one generated by a periodic potential, on this steady- 
state. We assume that the periodic potential has been 
switched on very slowly, so that in the absence of the 
initial quench, the system reaches the ground state in 
the presence of the periodic potential. The results are 
presented in the parameter regime where the periodic 
potential is (dangerously) "irrelevant" so that a pertur- 
bative RG approach remains valid at arbitrary length- 
scales. We find that infinitesimally weak potentials or 
mode-coupling can generate an effective-temperature and 
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cause the system to asymptotically thermalize. In addi- 
tion, and somewhat unexpectedly a dissipation is also 
generated. Thus we find that the effective low energy 
theory at long times after the quench is a quadratic the- 
ory of thermal bosons with a finite lifetime. This asymp- 
totic thermalization and dissipation occurs because the 
low frequency and momentum modes can transfer energy 
to the higher energy modes (which are gradually elimi- 
nated in the RG procedure). These high energy modes 
thus act as a bath, providing thermalization and dissipa- 
tion. The dissipation thus shares features with Landau 
damping where plasmons acquire a finite lifetime^ It 
also shares similarities with turbulent systems that are 
known to exhibit the well known Kolmogorov cascades 
where energy is passed down from large scale structures 
to small scale structures. 33 We believe that this type of 
behavior and mechanism for thermalization, unraveled 
in a controlled way on the particular system we study, is 
quite general. 

A short account of some of these results was given in 
Ref. lU In the present paper we give a detailed deriva- 
tion of the above results and in particular discuss the 
general derivation of the RG equations using the Keldysh 
approach which may be adapted to study other types of 
out of equilibrium bosonic systems. We build on the re- 
sults of Ref. HH to derive new results. In particular we 
show that the RG flow can be significantly modified by 
different quench protocols. We also examine the behav- 
ior close to the critical point where the external potential 
becomes relevant. In this limit various peculiarities such 
as singular behavior of the expansion in spatial and tem- 
poral gradients of the bosonic field are encountered. In 
addition we present a detailed discussion of the viola- 
tion of the fluctuation dissipation theorem of the post- 
quench system by studying the frequency dependence of 
two point correlation functions showing that even though 
the system appears "thermal" in low energy scales, the 
crossover from low to high energies is still complicated. 

It should be noted that some RG based approaches to 
study quenches already exist in the literature. For exam- 
ple the flow equation method was used to study quenches 
in a Fermi liquid^ and in the quantum sine-Gordon 
model, 34 where for both cases a long-lived prethermalized 
regime was found. RG was also used to study the a classi- 
cal two-dimensional sine-Gordon model after a quench^ 
Here the effective temperature due to the quench was 
found to generate a dynamical vortex binding-unbinding 
transition, but unlike the quantum ID problem studied 
in Ref. 3l| and this paper, no dissipation was generated. 
Furthermore, the system that we study also shares many 
common features with a Luttinger liquid subjected to a 
nonequilibrium noise source^ where the particular form 
of the noise in our system arises due to the out of equilib- 
rium occupation of the bosonic modes after the quench. 
It is thus interesting to make a connection between the 
quench problem where in the asymptotic the system is 
free, and systems where an external noise source is con- 
stantly imposed on the system, and is found to exhibit 



similar properties^ 

The paper is organized as follows. In Section|H]the ba- 
sic model and notation is introduced and the equilibrium 
properties of the relevant response and correlation func- 
tions are discussed. In Section HTT1 the interaction quench 
in the Luttinger liquid, and the properties of the resultant 
nonequilibrium steady-state is presented. In Section IIVI 
the periodic potential is introduced and the perturbative 
RG equations are derived. In Section [V] the solution of 
the RG equations are presented, and in Section PVTl the 
resultant low energy theory near the fixed point is dis- 
cussed. Finally in Section fVIII we present our conclusions 
and discuss open questions. 



II. MODEL 

The Hamiltonian for interacting bosons in a periodic 
potential is 



H = Hn 



Hn = — 



V, 



2 TV 



dx 



K (TrU(x)f + 1 (d x <fi(x)f 



v.. 



•— 2 / dx cos (7^(2;)) 
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(2) 
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where Hq is the quadratic part which describes the Lut- 
tinger liquid or long lived sound modes. The density 
of these modes is p — —d x (j)/w, whereas n = d x 9/ir 
is the variable canonically conjugate to (f>. V sg repre- 
sents the periodic potential whose most important ef- 
fect is a source of backscattering which can localize the 
density modes via the well known Berezenskii-Kosterlitz- 
Thouless (BKT) transition^ 

It is convenient to represent the fields 4>, 9 in terms of 
bosonic creation and annihilation operators (& p ,&|,),— 



4>{x) = -(N R + N L )- 



x l e -a\p\/2-ip* ( ft t 
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b- P ) 



ITT 

T 
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p#0 



2vr 
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(4) 



( X ) = {N R -N L )— + - 



E 



2ir 



1/2 



1 -a\p\/2-ipx 
\P\ 



where a 1 is an ultra-violet cutoff. Thus, 
H = ^u\p\r}ln P 

p^O 



where 



and e 2/3 



i 

K,u 



tip = cosh j3b p + sinh /3b 
- cosh /3b_ p + sinh j3b p . 

v F /K. 



(5) 



(0) 
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Since we are interested in nonequilibrium dynamics, we 
will use the Keldysh formalism^ where </>_/+ will denote 
fields that are (time/anti-time)-ordered on the Keldysh 
axis. Further, it will be convenient to define quantum 
(<j) q ) and classical fields (<fr c i), 



1 

7T 



cl T <t>q) 



(9) 



A. Correlation functions in equilibrium 



define the correlation function, 



C$ t {xt,yi) = -± (c^-(xt,yt') + Cjj+(U, v t')) (18) 

'] 



COS 



-ie 



K eq tan 



1 ( u(t — r ) + x — y 



+K eq tan 



l fu(t-t')-(x-y) 



(19) 



The two-point functions that are directly influenced by 
the periodic-potential, and therefore of interest to us are, 

Cff(xt,yt') = (e^iOM^-n-Mf.O) (10) 

In order to compute the above correlators in the absence 
of a periodic potential (g = 0), or perturbatively in g, 
some useful identities arc 

(VdMvtitfJ)) = Ve- <u|pl(t -° coth (11) 

(VclipMiP'^)) = S^eit-t^e-^-^ (12) 
(t? 9 (p,t)^(|/,0> - Spp>0V - t)e- iu ^-^ (13) 

where T is the temperature of the bosons. 
At T = and 5 = we find, 

C7r(xt,yt') = 



- 7 ^ zsign(£ — £') ^tan 



-1 u(t-t') + x-y | tan -l u(t-t')-(x-y) l 



(14) 



citify?) 

fin ~ — +ln 



■ | ,„ y •■ ■■■ ■■■■ ■' ■-' ,, ( a ,- w _„(t-t'))^ + ° 2 j 
^signtt-Oftan- 1 "C-''H*-i/ +tan -i (^(t-t')-(x-y)) ] 



(15) 



K \ ln V(x-y + u(t-t'))^+n2 | ^ v /( x - !< - 11 (t-t'))2 + 



+ln - 



-1 ».(t-t') + x-y 



+tan 



-1 ^(t-t')-(x-y) j 



(16) 



Thus all correlators exhibit the typical power-law decay 
of a Luttinger-liquid, with an exponent K eq where 



K eq = 



j 2 K 



(17) 



The oscillating factors in the above equations arise due 
to the Keldysh time-ordering and in equilibrium have 
the right structure so that the well-known fluctuation- 
dissipation theorem (FDT) is obeyed. To see this let us 



and the response function 

C^(xt,yt') = ^ (c^(xt,yt') - C^+(xt,yt')) (20) 
= -6(t-t') x 

_ Keg [l n \/(x-y + u(t-t'))*+<** | ln V(x-y- U (f-f'))^+c,^ j 



sin 



K eq tan 



_ 1 ( u(t — if) + x — y 



+K eq tan 



_! fu(t-t')-(x-y) 



(21) 



In Fourier space C(q,w) = f_ x dt J_ x dxe lu}t ~ lqx C(x,t) 
are in general complicated to compute. However we 
will in the subsequent sub-section only highlight how the 
FDT works, so that it is easy to follow how it is violated 
in the post-quench situation. 



B. Temperature from the Fluctuation Dissipation 
Theorem 



The FDT implies that 



C&(s,cj)= coth 21m [C&(q,uj\ 



T=0 



signH2Im [C&(q,uj\ 



(22) 
(23) 



Using ln(a + i(ut ± r)) = \n(y/a 2 + (ut ± r) 2 ) + 
itan -1 (placing the branch-cut of the logarithm on 
the negative real axis), we may write 



j>00 pOD 

C'^ l p(q 1 w) — — 2i / dt dr cos(wi) cos(qr) 
Jo J-oc 



-K ea In - 



-K ea \n- 



+ C.C. 



(24) 



/>00 />00 

2Im[C00(<7, u>)} — 2i / dt dr sin(cji) cos(qr) x 

JO J -oo 

!_ y Kcq , n ^±1) _ Keq ln ^M^l _ c c J (25) 
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Thus the FDT C K = 2sign(w)Im[C i? ] implies 

dr sin(oji) cos(qr) x 



sign(w) / dt 
Jo 

1 

2i 



K c „ In ■ 



+ i(ut + r) 



-Kealn- 



dt 



dr cos(wi) cos(gr) x 



1 r 

2 



o 

e Ae « ln s Ke<lln ° +c.c. 



= 

The above implies the following ought to be true, 



(26) 



Re 



dt 



10 



drcos(qr)e 8 l w l'x 



-K Ba In - 



= 



(27) 



By analytically continuing it —¥ t, it is straightforward 
to see that the expression in the square-brackets is purely 
imaginary thus proving the FDT at T = 0. 

In the next section, when we study the long time be- 
havior of the response and correlation after an interaction 
quench in the Luttinger liquid, we will find that the FDT 
is violated due to the appearance of a new nonequilib- 
rium exponent K neq which governs the power-law decay, 
while the oscillating factors are still associated with the 
equilibrium exponent K eq . 



III. INTERACTION QUENCH IN THE 
LUTTINGER LIQUID: PROPERTIES OF THE 
QUADRATIC THEORY 

Let us suppose that the system at time t < is a Lut- 
tinger liquid with interaction parameter Kq and velocity 
Mo, and therefore described by the Hamiltonian 



Hi 



2tt 



dx 



K {irIL(x)Y 



p^O 



(28) 
(29) 



We will consider the case where at t = there is an inter- 
action quench from Kq — > K so that the time evolution 
from t > is due to 



H f = — I dx 



2tt 



K(irTL(x)f + ±(d x <f>(x)) 2 



^u\p\lll P 

p^O 



(30) 



To preserve Galilean invariance (which is not necessary 
for the formalism), we assume u = vf/K,uq = Vf/Kq. 
Note that, 



,t 



cosh — sinh f3 
— sinh P cosh /3 

cosh /3q — sinh flo 
- sinh )3q cosh /3q 



Tp 

r-p 

Vp 



(31) 
(32) 



where e" 2/3 ° = K ,e~ 2/3 = K. The quench for K = 1 
was studied in Ref. [15|, and more general interaction 
quenches were studied in Refs. and HH. However, the 
distinction between response and correlation functions 
were only first identified in Ref. l3li Here we give more 
details of the results that appear in Ref. HH, and in ad- 
dition discuss the crossover behavior from low to high 
frequencies. 

Let us define the functions 

f(pt)= cos(u\p\t) cosh/3o 

-i sin(u|p|i) cosh(2/3 - f3 ) (33) 
g(pt)= cos(u\p\t) sinh 0q 

+i shx(u\p\t) sinh(2/3 - (3 ) (34) 

which determine the time-evolution after the quench [t > 
0), 

bl(t) + b. p (t)= (f*(pt) - g(pt))4(0) 

+ (f(pt)-9*(pt))v- P (0) (35) 

bl(t)-b- p (t)=(r(pt)+g(pt))r)l(0) 

-{f(pt)+g*{pt))n-p{Q) (36) 

Using the above, the basic expectation value for the cjy- 
fields after the quench can be easily worked out to give, 



4 y^oo ipi 



cos(p(x — y)) x 



\v\ 



— ^l + -^2jcosu|p|(t-i')coth — 
-iK smu\p\{t - 1')] 



(37) 



where T denotes the temperature of the Luttinger- 
liquid before the quench. Above, terms that oscillate 
as e -Hp|(*+* ) have been dropped. This is because such 
oscillating terms give an over-all decay to the correlators 
of interest defined in Eq. (|10p . Since we are ultimately 
interested in the long-time limit rather than the tran- 
sients, these terms are not important for us. Further, 
this approximation also leads to a significant simplifica- 
tion as one is now dealing with a nonequilibrium steady- 
state problem, defined by the following basic retarded 
and Keldysh Green's functions, 

G R (xt, yt')= -i6(t - t')([ct>(xt), <j>(yt')]) 

= -i{<f, cl (xt)cj> q {yt')) (38) 

G A (xt, yt')= i6(t' - t)([0(xt),<f>(yt')}) 

= -i{<jy q {xt)^ cl {yt')) (39) 

G K (xt,yt')= -«<{^(art),0(y*')}> 

= -i(Mxt)Myt')) (40) 

which in Fourier space acquire the following form at T — 
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G R (q,u) 
_ ttK 



■kKu 



(lu + i5y — u 2 q 2 
1 



1 



uj — u\q\ + i5 U) + u\q\ + iS 

.2 



.?r / , K 2 \ sign(w) 
G K (q,u>) = -JyA ( 1 + — I -, , x 



K 2 



[S(u> — u\q\) — S(uj + u\q\ 



2K 



1 



K 
K 



2 ; sign(w) [Gr - G A ] 



(41) 

(42) 
(43) 



Note that the retarded Green's function Gr depends only 
on the final Hamiltonian, and therefore is not sensitive 
to the quench. The Keldysh Green's function Gk on 
the other hand contains information about the occupa- 
tion probabilities of the bosonic modes, which can be 
far from thermal equilibrium due to the quench. Thus 
Gk depends on the properties of both the initial Hamil- 
tonian (via Kq) and the final Hamiltonian (via K). In 
equilibrium K = Kq, and the FDT (for T = 0) G K = 
sign(w) [Gr — Ga] is recovered. 



The power-law decay is determined by K neq which is 
a memory dependent exponent as it explicitly depends 
on the interaction parameter Kq before the quench. On 
the other hand, K eq depends only on the interaction 
parameter of the final Hamiltonian, and characterizes 
the equilibrium T = properties of the system. Fur- 
ther, K neq > K eq , so that the power-law decay is always 
somewhat faster in the nonequilibrium steady-state. The 
faster decay occurs both for (e 1 ^ 1 ) e -1 ^ 2 -') and the dual 
^s(i) e -t0{2)y j n sense, the effect of a quench is 
similar to a temperature, however the system remains in 
a critical state. In the next section, we will find that 
as a consequence of this, the periodic potential is always 
less relevant for the nonequilibrium steady-state prob- 
lem, with the critical point shifting to smaller values of 
K . Similar power-law decays with nonequilibrium ex- 
ponents can also arise in open systems subjected to a 
nonequilibrium noise source such as 1/f noised 

It is interesting to observe that two different quench 
protocols can lead to the same nonequilibrium steady- 
state, at least for a case where the steady-state is de- 
termined by the behavior of the correlators. To see 
this, for simplicity set 7 = 2 (so that K eq = K). Then, 



A. Correlation functions after the quench 

It is convenient to define a nonequilibrium exponent, 

(44) 



K neq = \K {1 + K 2 /K { 



2 lTr o) 



which we will show below represents the new power-law 
decay of the correlations in Eq. (fT0|). Interestingly both 
the equilibrium exponent K eq defined in Eq. (|17|) and the 
exponent K neq affect the unequal time correlators, which 
are found to be, 



ii^sign^-Oftan- 1 "('-0 + *-« +tal r' v(t-t> )-(*-y) 



(45) 



GU(xt,yt') 



-A', 



ejln 



^signCt-Oftan- 1 Mt-t') + *-v ^an' 1 ("('-'')-(*-«)) " 



(46) 



^K^ [tan- 1 " (f - f 2 +X '" +tan- 1 "('-''>-(*-»> ] 

(47) 

The above agrees with th e eq ual time (t — t') correlators 
for Kq = 1 studied in Ref lla. 



K„ 



K 



(K - Kq) 

2K a 



(48) 



From the above equation one may see that two different 
Kq may lead to the same K eq and K neq . For example 
K eq = l,K neq = 2 can be obtained for a quench from 
Kq —> K where K = 1, whereas the initial interaction 
parameter Kq can take two different values Kq — 2 ± \/3. 
This behavior simply reflects the fact that when K = Kq, 
the system being in equilibrium, K neq = K. On the other 
hand K neq increases with respect to K for both types 
of quenches, one where Kq > K and the other where 
Kq < K. For both these cases, the system is driven 
out of equilibrium, giving rise to a faster decay than in 
equilibrium. 

As before, we now discuss the FDT ratio for the 
two-point functions. (Note that often our convention will 
be to express length-scales in units of a and energy-scales 
in units of u/a). The Keldysh correlation function (de- 
fined in Eq. (fT8|) ) is found to be 



CW-M) = ~icos 

I i ut — r 

+K eq tan 



K eq tan 



l fut + r 



A„ 



a 2 + (ut — r) 2 



a 2 + (ut + r) 2 



A"„ 



(49) 



and the retarded correlation function (defined in 
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Eq. ([20])) is found to be, 
<7&(r,t) = -0(t)sin 



i4T e q tan 



--fT e( j tan 1 



ut — r 



a 2 + (lit — r) 5 



a 2 + (ui + rf 



(50) 



In general the Fourier transform of the above expres- 
sions for C K ' R need to be calculated numerically, and 
were briefly discussed in Ref. [U In the q — 0, uj — >• 
limit however, analytic expressions for can again be 
obtained. In particular 



C^(q = ) oj = 0) = -iI T o 
where (setting u = 1) 

^to = 
2a 2 



(51) 



2K».,-i(K neq -l)B 



K cq +K neq K neq -K e 



B(x, y) being the beta function. Similarly, we find 

21m [C^(q = 0,u -»■ 0)] = -iul n0 (53) 
where (setting u = 1) 
A70 = 



_2*-«- 1 (tf neg -l)B( 



2 K "-- 2 (^„ e9 -2)B 



2 ' 2 



2^»..-a(2f ne , - 2)B ( K -«+ K <"* , K ^-^ q -2 



(54) 



Note that in equilibrium K eq — K neq , and It,t)0 = 0. 



B. Violation of the quantum FDT and 
zero-frequency effective-temperature 

Even though the quantum FDT is not obeyed, one 
may define an effective-temperature T e ff.o in the low- 
frequency limit as follows, 



C&(? = 0,w = 0) 2T eff , 



21m CK(q = 0,u^0) 



(55) 



where we find that the effective-temperature (in dimen- 
sions of u/a) is 



T, 



Oil 



e//,0 



TO 



21, 



1/0 



2K 



(56) 



(::(/ 



As is typical of nonequilibrium systems, this effective- 
temperature depends on the correlation function being 
studied as it certainly does not characterize the low fre- 
quency properties of the simpler correlators in Eq. (|43p . 
Moreover, this temperature or equivalently the noise cor- 
relator has a complicated frequency dependence. How- 
ever, as we will show by doing RG, the low- frequency 
limit of the noise or the temperature has important phys- 
ical consequences, as it changes the long-time and dis- 
tance behavior of correlation functions by causing them 
to decay exponentially fast (rather than as a power-law 
with exponent K neq ). 

The crossover from uj <C ^e//,o to uj 3> 
T e ftfi is illustrated in Fig. Q] which plots the ratio 



2ImC«,(g=0,wy 

of uj coth 



2T„ 



comes uj coth 



In equilibrium, this ratio takes the value 
— which in the high-frequency limit be- 



2T, 



(52) frequency limit is uj coth 



|w|, and in the low- 



21: 



-> 2T, 



e//,0- 



The 



plot shows that the nonequilibrium system shows a slower 
crossover to \ui\ with increasing frequencies than the equi- 
librium system, indicating that the occupation of higher 
energy modes decays slower than exponential. It should 
be noted that similar noise with a complicated cross-over 
behavior from low to high frequencies was studied in open 
and driven systems near quantum critical points^r— 
The low-frequency limit of the noise was found to cut 
off the power-law decay of critical fluctuations, and to 

also cause a classical ordering-disordering phase transi- 
tionj 38 i 39 

Another measure of the violation of FDT is to extract 
the momentum dependence of the zero-frequency tem- 



perature, 2T ef f fi {q) 



2IraC? 



This quantity is 



, /,(<7.«->0)- 

plotted in Fig[5]and shows that the shorter the distances, 
the higher is the effective-temperature, unlike in equilib- 
rium where all length scales are associated with the same 
temperature. 

At this stage, a peculiarity of the zero-frequency ef- 
fective temperature, namely that it docs not vanish as 
K neq — > K eq should be noted. In this case, the first 
equality in Eq. (|56|) shows that the effective-temperature 
is a ratio of two quantities, both of which go to zero. 
However, the limit approaches a finite value. The origin 
of this result that the temperature approaches a finite 
value as the quench becomes smaller and smaller is due 
to the singular form of the equilibrium distribution func- 
tion coth(cj/2T), which probably persists even for the 
weakly nonequilibrium problem. The singular form of 
coth(w/2T) implies that the limits u> — » and T — > 
do not commute. For one order of limits the answer 
is divergent, and for the other it is 1. In defining the 
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zero-frequency effective-temperature, we have implicitly 
taken the frequency to zero first. This result also sig- 
nals that the energy scale T e ff$ cannot be used as a 
good measure of the energy stored in the system due to 
the quench. This energy is expected to be distributed 
in a rather complicated way among all frequency modes. 
However, by doing RG we will show that T e //,o acts like 
a regular temperature when studying two-point correla- 
tion functions as this energy scale causes the correlations 
to decay exponentially fast (as compared to a power-law) 
at long times after the quench. The size of the quench 
\Kq — K | on the other hand is inversely related to how 
long one has to wait to see this behavior. 

Besides an effective-temperature, perhaps a more sur- 
prising result, is a generation of friction. This effect ap- 
pears at this stage as a non-zero slope of 2Im[C^] a 
—iijco. In the next section when we do RG, we will 
show that 77 corresponds to a finite-lifetime of the bosonic 
modes, an effect which is distinct from the generation of 
an effective-temperature. Thus we will find that even 
though the system we study is closed, and in equilib- 
rium is characterized by long-lived bosonic modes with 
•q = 0, the quench together with the mode-coupling aris- 
ing due to the periodic potential gives rise to additional 
scattering which generates an 77. The finite lifetime of 
low-frequency bosonic modes implies that there is a flow 
of energy from low energy scales to high energy scales. An 
alternate, but simpler example of this phenomena is the 
decay of collective modes of a system of one-dimensional 
weakly interacting fermions via the creation of particle- 
hole excitations, where the fermions are in a nonequilib- 
rium state due to an initial quench^ 



IV. DERIVATION OF THE RG EQUATIONS 

In order to derive the RG equations, it is convenient 
to write the Keldysh action for the steady-state, 



J K 



= jT>[<j> c 



J(S +S sg ) 



(57) 



where So is the quadratic part which describes the 
physics at long times after the interaction quench, 





G r> 



G A 

— 1 { ' f" — 1 



where, 



t) (£■ O) 



j Gk Gr 
\G A 



(58) 



(59) 



Using the expressions for Gr^a,k in Eqns. (HTj) and (|42l) . 
Eq. (|58[) may be written in the following manner, 



So 



_L » 

ttKu y{uJ + iS) 2 -u 2 q 2 

4> q (q,u) 



{lu — iS) 2 — u 2 q 2 



(60) 



where 5 = 0+, and the modes are long-lived. Note that 
the above action though obtained for the specific case 
of an interaction quench in a Luttinger liquid, has the 
same generic form as any Luttinger liquid subjected to 
a nonequilibrium noise source. The details of the noise 
determines the particular form of the coefficient of the 
(f> 2 -teim in the action<&2& In fact nonequilibrium noise 
can also arise in open systems driven out of equilibrium 
by current flowj^ 

We now discuss how the above free theory is affected 
by a mode-coupling term due to a periodic potential. The 
action corresponding to this is 



S a 



gu 



dx J dt [cos — cos 70+] (61) 
We split the fields into slow and fast components, 



4>cl,q(xt) = (f>5jxt) + (/>>, (Xt) 



(62) 
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and integrate out the fast components. We will out- 
line two different procedures for integrating out the fast 
modes. In one we impose a hard cutoff, where the fast 
fields arc defined as, 



poo 


duj r A ' 


dq 


J — oo 




2n 


did 


f 


dq 




'A>|g|>A' 


2^ 



^iqx—iujt 



</>q tC i (q, wj)64) 



and A/ A' = e dln(l \ The results in Ref.|U 

were presented 
far from the critical point using the above RG scheme. 

The second way of integrating out fast modes is a 
scheme outlined by Nozieres and Gallet,— where the slow 
and fast correlators are defined as follows, 



Go,. 



^A-dA.A 



(65) 



Thus the fast correlator may be obtained from taking a 
derivative of the slow or full correlator, 



r> -H.\ d ° A 
"A— <2A,A — ~JX~ 



(66) 



where A = a -1 . Note that when doing RG in real time, 
all cutoffs should be imposed only in momentum space, 
as imposing cutoffs in time and varying them during the 
RG flow lead to inconsistencies such as a violation of 
causality. 

The Nozieres-Gallet scheme is a more consistent way to 
deal with the cutoff, and finally the results in this paper 
will be presented using this method. However even the 
hard cutoff scheme outlined in Eqns. (|63)l and (|64l) and 
used in Ref. |3f| gives qualitatively similar results. Fur- 
thermore, in this paper we will use the Nozieres-Gallet 
scheme to study the behavior close to the critical point. 

Now our task is to expand Zk in powers of g, inte- 
grate over the fast modes, and rescale the cutoff back to 
its original value (dx(dt) ~ > -p-dx(dt)). Up to quadratic 
order in g we find, 

Z X = fv[fo*<\ e iS ? 

Xe -T<(^) 2 > (cos(7^<(l)) - cos(7^(l))) (67) 



, , -9U ( A 



dxdt 



2 2 

9 u 



J dxidtidx2dt20(ti — t%) 



(cos ( 7 (0<(1) + 0>(1))) cos ( 7 (0<(2) + 0>(2))))> 
— J dxidtidx2dt 2 0(t 2 - h) 

(cos ( 7 (^(1) + ^(1))) cos (7(^(2) + ^(2))))> 
a 2 u 2 r 

/ dXldtldX2<lt2 [ fl (*2 ~ *l) + - *2)] 

(cos ( 7 (0<(1) + </>>(!))) cos ( 7 (0<(2) + 0>(2))))>] 



The above may be re-exponentiated using the cumulative 
expansion (e v ) ~ e^ v ' + ^^ v to obtain 



Zk 



(69) 



The expression for SS may be simplified by dropping 
terms that are proportional to e i ^ ± ^ +i ^ ± ^ as they are 
more irrelevant than terms such as e l ^ ± ^~ l ^ >± ^ 2 \ Fur- 
thermore, using the fact that 



(sin [ 7 0>(I)-70>(2)]> =0 



(70) 



all the terms containing sin- functions vanish. Moreover, 
using cos(/> =: cose/) : '' 2 we find, 

mi / A \ 2 t t ^ 2 ^A^, 

5S ■ 

a' \ A' J J J 

(71) 



a 2 \M J J J 

x (cos(70<(l))-cos( 7 ^(I))) 



+i 



2 2 

g u 
2a 4 



dx\ J dt\ J dx2 J dtz6(t\—t2) 



: cos( 7 (0<(l) - 0<(2))) : e -¥<W-UW-(2)) 2 > 

i - e -4 ((4>> l ) 2 )+4 <w>(i)-0>(2)) 2 > 

g 2 u 2 f 

/ dxidtidx 2 dt20(t2 - h) 

: cos (7(^(1) - ^(2))) : <W+(D-<M2)) 2 > 
1 - e -4 ((<t>>if)+4 <(0?(i)-0?(2)) 2 ) 

J dx 1 dt 1 dx 2 dt 2 [6{h - h) + 0(h - h)] 
cos( 7 (0f(l)-<£<(2))) :e -^<W-(i)^+(2)) 2 > 

1 - e -4 <W c > 1 ) 2 > + ^<(0>(!)-^(2)) 2 



(72) 



(73) 



2 2 

.9 u 

l t 

2a 4 



(74) 



where ((<f>± (I) — 4>± (2)) ) involves averaging all (both 
slow and fast) modes. 
Eq. fl7T} implies, 



m = .9(A) 



A 



e 4 



(75) 



g(A) (-) exp 



-L-K (l + K 2 /KZ)HA/A') 



(76) 



Note that 

1 - e -T<Wc > ! ) 2 > + T(Wa(l)-^ > (2)) 2 ) 



(77) 



>J Next we introduce the center of mass R = Xl + X2 , T m = 

(68) tl t t2 and relative coordinates r = xi — x 2 ,r = ti — t 2 , 
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and expand the 0(g 2 ) terms in powers of r = X\ — X2, t 
ii - 1 2 - Using, 



:cos( 7 (0<(l)-0<(2))) :~ 

1 " \ [(rd R + rd Tm ) <f> cl + (rd R + rd Tm ) 9 ] 2 (78) 
:cos( 7 (0<(l)-0<(2))) :~ 

1 " \ [(rd R + rd Tm ) <f> cl - (rd R + rd Tm ) 9 ] 2 (79) 
:cos( 7 (0<(l)-0<(2))) ^ 

1 - -j l(rd R + rd Tm ) <\> d + 2</> q ] 2 (80) 



we note that the terms with 1 and purely classical fields 
cancel. We regroup the remaining terms and find the 
following corrections to the quadratic part of the action, 
(the correction to the cosine term is already given above 
in Eq. (75}) 



x*)/(2i)) 



It = 



dr / dt 



Re 



Im 



e -i-<(<Mt,r)-<M0,0)) 2 ) pa/ng 

dr / dtt 



— oo <J — CO 



- J T ((<t>-(t,r)-4> + (0,0yf) pa,ng 



I K = I dr dt (r l - t z ) 

l-oo JO 



Im 

In : 

Im 



-X_((0_(t,r)-0 + (O,O)) 2 ) pa/ng 



dr / dt(r z +r) 

oo JO 



e - J 2-((<P-(t,r)-<P + (0,0)) 2 ) pa/ng 



(87) 



(89) 



(90) 



where F a arises due to the hard cutoff scheme and F ng 
is due to the Nozieres-Gallet scheme. In particular, 



S Q = JdRj d(uT m )-L [4, q (d 2 R ~ d 2 uTii _ 



5u 



+0d (pl - d 2 uT J <i> q + -4> q {pi + d 2 uT J fa 



AT effV K ( K 2 
+l ^ 2 ~2K\ l+ K 2 



The above implies the following RG equations, 



(81) 



e -i-<(0-(^t)-0+(o,o)) 2 > = 

& i [K cq tan" 1 ( i±f ) + K eq tan" 1 ( i=t )] 



(91) 



* / K, 



eq 



2 \K„ 



n eq 



1 {sin(A(t + r)) + sin(A(< - r))} 



(92) 



a 2 + (t + r) 2 a 2 + (i-r) 2 
a(i + r) a(t — r) 



neq 



a 2 + (t + r) 2 a 2 + (t-r) 2 



(93) 



dg 

dK- 1 
din 2 



2-^/f (l + X 2 /^ 2 ) 



4a 4 V 2 



Kudlnl 4a 4 V 2 



/Y 



din/ 



= V 



ng 2 Ku ( 7 2 \ 2 Kq 



K 
K, 

d(T ef fV) OT _ , tt 5 2 U 2 ^ 2 ^ 2 



din/ 



2a 4 
2T e// ?7 



2 j "T 1 



(82) 
(83) 
(84) 

71^ (85) 



4a 4 



y) It 



(86) 



where (defining Re(x) = (x + x*)/2,Im(x) — [x 



V. RESULTS IN THE GAPLESS PHASE 
A. General structure of equations 

In this section we study the consequence of 
Eqns. ([82]), (53}, dHJ), (85} and (86). Eq. (82} shows 

2 

that there is a critical point located at K neq = ^-Kq(1 + 
K 2 /Kq ) — 2, and therefore at a value of K which is dif- 
ferent from the equilibrium critical point at K = 8/ 7 2 . 
Since K neq > K eq , this critical point is always located 
at a smaller value of K , which is another way of saying 
that the cosine potential for the post-quench case is more 
irrelevant. 

Eq. (53} represents the flow of the interaction constant, 
however depending on the quench protocol (and hence 
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the values of K eq and K neq ), the flow can be significantly 
different from the equilibrium flow. For example, while 
Ik > in equilibrium, implying a decrease of the interac- 
tion parameter in the presence of the periodic potential, 
for the quench problem Ik can become negative, and 
also diverge at the critical point. This will be discussed 
in more detail below. Eq. (|84|) is the flow of the velocity 
which occurs even in equilibrium, and is primarily due to 
the cutoff procedure employed here which does not pre- 
serve Lorentz invariance. The effects of this are small, 
and in what follows, we will ignore it. 



Eq. (|85|) shows that a dissipation or a finite lifetime 
of the bosons r\ is generated, changing the low-frequency 
properties of the bosonic system qualitatively. Further, 
I n (which is the rate at which rj increases with flow) di- 
verges at the critical point implying a diverging dissipa- 
tion. The origin of this divergence is similar to the di- 
vergence of Ik briefly mentioned above. This divergence 
can have two possible causes: i) the flow is derived using 
the original (athermal state) correlation functions. Since 
a finite temperature and dissipation is generated, this 
might regularize the divergence close to critical point; ii) 
the renormalization of the coefficients rj and K results 
from a gradient expansion of the second order results. 
Such a divergence in the correction might indicate that 
this expansion breaks down. This could mean a non- 
analytical behavior at low energy. Unraveling this point 
is a challenging question. 



Finally Eq. ([86j) shows that a constant term is gener- 
ated for the strength of the noise correlator in the zero- 
frequency limit, which can be interpreted as a product 
of dissipation and effective-temperature. In general this 
noise is expected to have a complicated frequency de- 
pendence as was highlighted for the g = case in the 
previous section. To determine how this frequency de- 
pendence evolves with RG is a daunting task, and in this 
paper we will only discuss the effects of g on the low- 
frequency part of the noise spectrum. 

In both the hard cutoff scheme as well as the Nozieres- 
Gallet scheme, the important result that It^ = when 
K eq — K neq is recovered. Further both schemes reveal 
the peculiarity of diverging Ik and I v at the critical point 
Kneq = 2 for the nonequilibrium problem. It should be 
noted that in deriving the above RG equations we have 
used the correlators for T e ff = rj = 0. The more consis- 
tent way to treat the problem is to evaluate the corre- 
lators for non-zero T e ff and t). However this makes the 
problem quite difficult, all the more so because the full 
frequency dependence of T e ff is needed. We argue that 
taking these effects into account only has a minor influ- 
ence on the RG flows far from the critical point. However 
the divergences in I n ,K indicates that near the critical 
point, a more consistent computation may be needed. 
This is clearly something left for future studies. 

In the Nozieres-Gallet scheme, the expressions Ix,n are 
found to be related in a rather simple way to Ir.nO i.e., to 
the q — 0, oj — » limits of the C K and lmC R evaluated 



for the quadratic theory in Section IlIIl In particular, 



K neq (K neq 1) 



K — 1 

~n — 1 

- 1 »- in n 



TO 



(94) 



I v = 8a 6 



VK.„ K nca -K c 



i eq 



2 ' 2 
(K^ - 3/2) 

-K ne q 3/2 



K n 



eq 



(95) 



where It v o are given in Eqns. (1521) and (|54|) . The ratio 
of these quantities leads to an effective-temperature, 



ng _ air _ (K neq - l)(K neq 
eff ~ 21, 



2K eq (K neq ~ 3/2) 



K, 



neq 



l 



K ne q 3/2 



T, 



e//,0 



(96) 



(97) 



with T e ff t Q being the non-interacting expression for the 
effective-temperature. Note that when K neq ^> 1, the 
cosine potential is more irrelevant, and the above expres- 
sions approach the non-interacting values. The above 
analytic expressions also show that I n is divergent at the 
critical point as l/(K neq — 2). 

The expression for Ik is more complex and given by, 



/ dx tan xcos(K eq x) (cos x ] 
Jo 

/ dytanysin^^yXcosy)^""- 2 
Jo 



K„ 



ir/2 



K„ 



dx tan x cos(K eq x) (cos x) 



x / dytanysuY(K eq y)(cosy) Kn "> 



K, 



eq 



Kn 



tt/2 



dx tan x cos(if e<? a;)(cos x) Kn "' 3 



l ne<j / Jo 

x / dy(tany) 2 cos(K eq y)(cosy) Kn,!q 



(98) 



(99) 



(100) 



K, 



eq 



K„ 



71-/2 



dx(tanx) 2 sm(K eq x)(cos x) Kn ° q 



'•-neq/ JO 

x / dy tan ysm(K eq y)(cosy) Kn ' !q ^' 2 



(101) 



Ik shows complicated structure, and in particular de- 
pending on the quench protocol can be positive, nega- 
tive or zero. For example, lK(K eq = 2,K neq = 7) — 
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FIG. 3. Change of sign of Ik at a certain value of K neq . Here 



— a 4 0.0022 showing that for this quench protocol, the 
periodic potential increases the interaction parameter, 
making the system more delocalized as compared to the 
ground state of a Luttinger liquid with interaction pa- 
rameter K. A plot of Ik for K eq = 2 and different K neq 
is shown in Fig [3] where a change of sign of Ik is found. 

Further, Ik generically diverges at the critical point 
K neq = 2 for all values of K eq except K eq — 1 and 
K, 



eq 



2. This divergence is seen from observing 
that when K neq = 2, the Eqns. (|M|) and (flUO")) above 
are in general divergent because tan(7r/2) = oo and 
(cos(x)) Kn "'~' 2 = 1. However these divergences can- 



cel in equilibrium when K e 



Kn 



2, and we find 



lK(K eq = K neq — 2) — a 4 (87r/32). In contrast at the 
special point K eq = l,K neq = 2, lK(K eq = l,K neq = 
2) = 0. 

In the subsequent sections we will study the RG flows 
in the gapless phase for two cases, one where the system 
is far from the critical point, and the other where it is 
close to it. For the latter we will highlight the effect of 
the above divergences. 



B. Flow far from the critical point 

The RG equations can be solved numerically. The flow 
of the various quantities for a quench from Kq — 3 to 
K = 5 and for 7 = 2 are shown in Figs. HI O H] and [7J 
For these values Ik remains positive, so that K decreases 
during the flow (Fig. [5]) . Since g is (dangerously) irrel- 
evant it flows to zero (Fig. 2]) , however for the nonequi- 
librium problem, g generates new terms such as a dissi- 
pation (Fig. [6]) and a zero-frequency component of the 
noise (Fig. [7]). The latter may be identified from the zero 
frequency limit of the following term in equilibrium, 



T]u> coth(oj/2T) 2r]T 



(102) 



as dissipation x effective — temperature. Note that since 
we have, from the RG, a quadratic action for the fixed 
point, this allows us to unambiguously define a tem- 
perature for the low energy modes. Indeed with such 



a quadratic action the fluctuation-dissipation theorem 
would be obeyed for all correlation functions involving 
the field 4> with the same temperature T as defined 
by (|102p . One thus sees that the RG procedure indi- 
cates that the system does thermalize. It is important 
to remember however that this statement only concerns 
the low energy modes. The frequency dependence of 
the Keldysh term is quite different from the standard 
r?wcoth(u;/2T). This means that as the frequency is in- 
creased there will be a complicated crossover between 
this thermal state at low energy and the athermal orig- 
inal state that controls the high energy behavior of the 
system. Another interesting question is what will hap- 
pen for the correlation functions of the dual field 9. The 
most naive expectation would be that such correlations 
are also controlled by the same temperature T. However 
since in the above action 9 has been integrated out, this 
would need explicit calculations of the 9 — 9 correlations, 
a quite complicated calculation. It will thus be interest- 
ing to see whether this is indeed the case or whether the 
dual field can have a different behavior. 

The flow of the effective-temperature can be studied 
most easily by substituting Eq. (|55"j) in Eq. We 
obtain, (after defining dimensionless variables T e ff — ¥ 
aT ef f/u,rj -> arj/u) 



dT, 



eff 



d\nl 



T, 



IS 



1 - 2T, 
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7T5 
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Ko 
2K 



K 2 \ I, 



K 2 
^0 



Note that the physical temperature T e fft 
differential equation 

d (T eff e 



I T 

-In l 



(103) 
obeys the 



lnl\ 



d\nl 



1 - 2T, 



eff 
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2K 



4 

K 2 



1 



K 2 
^0 



It 



(104) 



The above shows that the steady-state solution corre- 
sponds to a temperature, 



1 eff ~ 



K 



neq 



1 



K* — 2 
2K* 



(105) 



Note that T 



eff 



K* — 3/2 

neq *-V / \ "-''■neq 

differs from the value of It/(2I v ) in 



Eq. (|97p only because of the way we have defined it. The 
energy-scale which determines the decay of the correla- 
tions is the prefactor of the <fi 2 term in the action Eq. (|81|) , 
and is the combination T*^K^ eq / K* q . 

The renormalized dissipation e~ ln '77 for a quench from 
Ko = 3 to K > Ko is shown in Fig. [51 The behavior is 
non-monotonic as the size of the quench becomes larger 
and larger. The reason for this is that when K = Ko, r\ = 
0. At the same time when K ^> Ko, the cosine potential 
being more irrelevant, decays faster to zero, so that the 
renormalized r\ is also smaller for larger K. These two 
behaviors for \K — K \ *C 1 and K Kq imply a maxima 
in between. 
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FIG. 4. Flow of g with x = In/ for K = 3, if = 5, g = 0.1 FIG. 7. Flow of e~ 2x T eff rj with x = lnZ for K = 3, K = 5, 
and 7 = 2. g = 0.1 and 7 = 2. 




The generation of a temperature and dissipation in 
the low-energy theory via the RG procedure is shown 
schematically in Fig. [SJ For the out of equilibrium sys- 
tem, there is a gradual flow of energy from the low-energy 
(long wavelength) degrees of freedom and the high-energy 
(short wavelength) degrees of freedom. Thus when the 
latter are integrated out in the RG procedure, this flow of 




0.2 0.4 0.6 0.8 1 



x 

FIG. 6. Flow of e- x r 1 (x) with x = In I for K = 3, K = 5, 
g = 0.1 and 7 = 2. 



energy appears as a dissipation in the low-energy sector. 
Such a dissipation is also accompanied by a temperature 
such that a low-frequency classical FDT is obeyed. The 
above mechanism for thcrmalization arising due to an 
exchange of energy between long wavelength and short 
wavelength modes appears to be rather generic, and may 
even be recovered in a quench involving fermions^ In 
particular in Ref. HH the effect of weak interactions on a 
system of one-dimensional fermions that are in a nonequi- 
librium state due to an initial quench was studied using 
the random phase approximation (RPA) . The RPA anal- 
ysis revealed that the highly nonequilibrium fermion dis- 
tribution generated by the quench results in an enhanced 
particle-hole continuum. Thus for attractive interactions 
between fermions the collective modes were found to lie 
within this continuum, and were therefore found to be 
overdamped. 

It would be interesting to study other models with non- 
linearities to see if a similar generic fixed point is reached, 
and also how the temperature might depend on the type 
of non-linearities that exist in the Hamiltonian. Among 
the various non-linearities that one could think of, there 
exists in particular the band curvature that is inherent 
to a realistic one dimensional system. Such a band cur- 
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FIG. 10. Strength of the dissipation r\ for g — 0.05, Ko = 1 
and 7 = 2 so that K neq = 1(1 + K 2 ). The critical point is 
located at K — y/3. 



FIG. 9. A schematic of the RG procedure. The resultant dis- 
sipation and temperature generated in the low-energy theory 
(box) arises due to an exchange of energy with the high-energy 
degrees of freedom that act as a bath. 



vature leads in particular to (V0) 3 terms. Whether such 
terms could lead to similar effects is an interesting ques- 
tion. Note that in equilibrium, because of conservation 
laws, the cosine terms are more efficient in mixing certain 
modes in contrast to other non-linearities. In particular 
the curvature terms were found not to be able to fully re- 
lax a current in Ref l43l Whether similar differences also 
exist out of equilibrium is clearly a challenging question. 

Independently of these subtle points we think that the 
mechanism that we could obtain in a controlled way for 
this particular model using the RG procedure is quite 
generic. Fig. |H] indicates that provided enough mode cou- 
pling exists in a system, the system itself can act as a 
reservoir and bath for the low energy part of the degrees 
of freedom. This sub-part will thus acquire a classical 
behavior in the sense that it will get a finite dissipation 
and a finite temperature. It is important to contrast this 
with a standard equilibrium quantum system for which 
the temperature is the same irrespectively of the energy 
of the mode considered. Here the temperature can only 
be defined if the limit u — > is taken. The frequency de- 
pendence will in general be complicated and correspond 
to the crossover between the athermal and the (low en- 
ergy) thermal state. 



C. Flow in the vicinity of the critical point 

In this section we study the RG equations near the 
critical point and in particular highlight the effects of di- 
verging Ik.t] ■ Fig. [10] shows the renormalized dissipation 
for 7 = 2 and for a quench from Kg = 1 and a K chosen 
to be K > y/3. Thus these quenches include those that 
go across the equilibrium critical point at K = 2, and are 



always on the gapless side of the new nonequilibrium crit- 
ical point located at K = v3. Fig. HOl shows that as K is 
varied such that one approaches this new critical point, 
the dissipation diverges due to diverging I n . This could 
either signal a breakdown in the gradient expansion for 
the nonequilibrium problem and/or a drawback of our 
approximation of setting T e f f = r\ = in the two-point 
functions used in the evaluation of the RG equations. 

Neglecting the effects of temperature and dissipation, 
it is interesting to study the BKT flow near the quantum 
critical point. It is of course understood that the flow 
would be eventually cutoff by the temperature and/or 
the dissipation. Note that close to the quantum critical 
point K neq = 2, the temperature (|105[) becomes paramet- 
rically small while r\ diverges, so that the product rjT e f f is 
well defined. Yet, the parametrically small temperature, 
probably signals that as one approaches the Mott-phase, 
the low energy theory is no longer described by gapless 
thermalized modes. This is an important open question 
which we do not address in this paper. 

Near the critical point we may expand the RG equa- 
tions for g and K (Eqns. (|82p and (|83[) ) about K neq — 
2 + e. For concreteness let us set K eq = 3/2. Then the 
complicated expression for Ijc in Eqns. (j9"5)l . (j9"9")l . (|100[) 
and (jlOip reduces to 

I K {K eq = 3/2, K neq = 2 + e) = -a A [-0.54 



+ — / ax tan x cos x) sm — 
7 Jo V 2/ 

x f cos — J (3 + 3 cos x + cos 2x) 
' (cosx) c ^sin — ^ 



7T/2 



dx tan : 



3.r 



I cos — ) (2 + 9 cos x + 3 cos 2x) 



(106) 



where the first smooth numerical term comes from 
Eqns, (r9"8"|) and (|101l) . while the rest are terms that di- 
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e 



FIG. 11. Flow near the critical point K neq = 2 + e, neglect- FIG. 12. A plot of rj/g 2 vs T 2 where T = T* ff K* eq /K* q and 

ing the effects of temperature and dissipation. Note that the for quenches where 7 = 2, Ko = 3 and K > 3. 
boundary separating the regions in which the cosine is rele- 
vant or not is quite different from the equilibrium BKT one. 



verge as e — > 0. Thus for e <C 1, we may approximate 

Ik ^ - (107) 



where c is a positive number. Thus after a redefinition 
of g, the flow of g and e are found to be 



dg 
d\nl 
de 
dial Z 



„2 



(108) 
(109) 



The above imply that the flow equations are along the 
following lines 



9 2 A 
= A 



(110) 



where A is a constant. Some examples for different 
choices of A are shown in Fig. [TT] and should be con- 
trasted with the equilibrium BKT flow which are along 
the lines g 2 — e 2 = constant. 

Along the separatrix A = 0, the solution of the RG 
equations give 



e(l) 



1 + §e lnZ 



1 



3/2 



(111) 

(112) 



_l + (2|o) 2/3 lnZ_ 
(denoting eoj <?o as the initial bare values). 



VI. PROPERTIES OF THE FIXED POINT 
ACTION 

In this section we discuss the general properties of the 
steady-state resulting from the RG. The renormalized ac- 
tion is a quadratic theory of thermal bosons with a finite 



lifetime. In the low- frequency limit, the effective action 
is given by 



So= ^2 (^*/(9> w ) ^qil^)) 



_L 

txKu \uj 2 + irjuj 



UJ — ITjUJ - 



2 2 

u z q z 



2 2 

u q 



(113) 



Thus the 



correlators are, 



uKu 



uj' z — u 2 q 2 + irjU! 



x (2T e// r?) 



2K 
uK 



rfuj 2 + (uj 2 — u 2 q 2 ) 2 



(114) 



(115) 



The above correlators should be contrasted with those 
of the nonequilibrium Luttinger liquid discussed in sec- 
tion III. While both effective theories are quadratic, the 
combined effect of a quench and the periodic potential 
is to give rise to inelastic scattering that broadens the 
bosonic modes by an amount given by rj. For the Lut- 
tinger liquid on the other hand (section III), the bosonic 
modes are long lived, but are characterized by a nonequi- 
librium occupation probability. Note that a finite lifetime 
r\ may also be generated for interacting bosons in a peri- 
odic potential which is in equilibrium but at a non-zero 
temperature, provided there are no special conservation 
laws4M5 

In what follows we neglect the lu 2 term relative to 
the r)uj term as we are primarily interested in the long- 
distance and long-time limits. Then the retarded corre- 
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lator is found to be, 

-1(^(1)^(2)) 
dq f doj 



2tt / 2tt 



= -0(h-hY- 



iq(xi-x 2 )-iu(t 1 -t 2 )QRfq^ ^ 



K 



2 - 1 2 ) 

Similarly the advanced correlator is, 

-i<^(l)&,(2)> 
dq f du 



e 4 " 2 i*i- f 2i (116) 



2tt / 2tt 



= -e(*2-tl)- 



g(xi-x 3 )-<u)(ti-t2)g!A^^ a; \ 



A" 



2 Vvih - h) 

while the Keldysh correlator is found to be 

-i [(^(1)^(2)} - <^i>] 



e ^| tl -t 2 i (n7) 



d<7 /" dw 
2tt 7 2tt 
A 



D iq(x 1 -x 2 )-iu)(ti-t 2 ) _ ^ 



A 



M) 



It, -t 2 l ar 



n 



X\ - X 2 



Erf 



y/rjKx! - x 2 )\ 

2uy/\h-t 2 \ 



(118) 



For ti — t 2 , the above reduces to 
-i [(<t>ci(x 1 ,t)<p c i(x 2 ,t)) - {(j> 2 cl )] 



2K 



1 + 



El 

K 2 



m ff K) 



(119) 

Eq. (|1 19|) shows that the equal-time two-point function 
C¥a decays exponentially in position, with a decay rate 
given by the effective-temperature, T e ffK neq /K eq 

(120) 



( e i<Pc(x) e - 



i<Pcl(v)\ 



An interesting question concerns the relation between 
the generated dissipation and the effective-temperature. 
For a fermi-liquid for example r\ ~ g 2 T '. It is interesting 
to explore to what extent our system mimics such a be- 
havior. Fig. [T2l shows the plot of r]/g 2 as a function of 
T 2 where T = T* ff K* (tq /K* q , the latter being the ap- 
propriate energy scale that determines the decay of the 
two-point functions. The * indicates that the values at 
the fixed point have been taken. Fig. [12] corresponds to 
7 = 2 and a quench from Kq = 3 to K > Kq. No- 
tice the coincidence of the plots for two different values 
of g indicating that the dissipation scales as g 2 . How- 
ever, for small quenches, the behavior is not fermi-liquid 
like as the dissipation increases as T@ with j3 > 2. For 
larger quenches, the relation between rj and T becomes 
non-monotonic (not shown) as T always increases with 
AT,whereas 77 eventually decreases as g becomes more and 
more irrelevant. 



VII. CONCLUSIONS AND OUTLOOK 

In this paper we have explored how out of equilibrium 
quantum systems can thermalize in the presence of mode- 
coupling terms. We have analyzed in detail a situation 
in which a system of interacting bosons is set out of equi- 
librium by a sudden interaction quench. In Ref. [3l| and 
the current paper we have shown, using a controlled RG 
procedure, how mode-coupling or interactions which give 
rise to nontrivial scattering between modes affect this 
nonequilibrium state. 

The main result is that even when the mode-coupling 
is "irrelevant" , it generates an effective-temperature and 
a dissipation. The generation of a dissipation indicates 
that thermalization eventually sets in by the exchange 
of energy between low energy modes and high energy 
modes. The flow of energy across different length scales 
is also found in nonequilibrium classical systems such 
as turbulence^ and the classic Fermi-Pasta-Ulam prob- 
lem^ Such a flow of energy can lead to cascades and uni- 
versal power-law behavior in the distribution function, a 
result which has also been recovered in recent studies in- 
volving nonequilibrium Bose-Einstein condensates^ It 
is clear that exploring further the connection between 
these systems should prove to be a very fruitful line of 
study. This also points to quite a general mechanism for 
thermalization in which mode coupling allows the sys- 
tem itself to act as a reservoir for the low energy part of 
the degrees of freedom (see Fig. [9J . This subpart thus 
acquires a classical behavior in the sense that it is char- 
acterized by a finite dissipation and a finite temperature. 
It is important to contrast this to a standard equilibrium 
quantum system for which the temperature is the same 
irrespective of the energy of the mode considered, here 
the temperature can only be defined if the limit uj — > 
is taken. The frequency dependence will in general be 
complicated and will correspond to a crossover between 
the athermal and the (low energy) thermal state. 

It is also interesting to compare the approach we have 
employed with more traditional methods to study dy- 
namics. When studying fermions, a natural route is to 
write a kinetic equation that systematically takes into 
account two and if needed three particle scattering pro- 
cesses. However for the Luttinger liquid which has a lin- 
earized spectrum, a description in terms of a kinetic the- 
ory fails. In particular a naive perturbation theory about 
the Luttinger liquid fixed point leads to divergent results. 
This has been attempted in the past (see Ref. |48| for a 
discussion on this point). While our approach is pertur- 
bative in the periodic potential, since the basic unit of our 
perturbation theory involves the correlators ( e M e - w< A) j 
we have explictly taken into account multiple scattering 
between bosons. 

There are of course many open problems. It is in par- 
ticular important to extend the analysis of the present 
paper close to the critical point and also to the region in 
which the mode coupling term is perturbatively relevant. 
Although we could obtain in the present paper some re- 
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suits close to the critical point, the whole RG procedure 
has to be made fully consistent by taking into account the 
temperature and dissipation to obtain a more complete 
description. Treating these effects in a self-consistent way 
is hard. Once such a scheme is developed, it would be 
interesting to study the strong-coupling part of the phase 
diagram. It is clear that numerical studies of these ques- 
tions would also be extremely helpful in this regime. 

Given the close analogies between fermions and bosons 
in ID, exploring the above physics in fermionic systems 
is an important direction. A first step was undertaken 
in Ref. [32] where an RPA study of a system of fermions 
that are out of equilibrium due to an initial quench was 
done, and overdamped collective modes were recovered. 
However, a more complete study is needed that takes 
into account backscattering interactions (which are not 
included in the RPA). Further, we have assumed that 
the term that generates mode-coupling was switched on 
very slowly. How the physics is affected by the rate at 
which the mode-coupling is turned on is also interesting 
to explore. 

Finally checking the generality of this mechanism, by 
introducing other non-linear terms, such as those aris- 
ing due to band- curvature, is of course a very interest- 
ing question. This is also directly relevant for a test of 
the present mechanism either in numerical simulations 



or in experiments. On a practical point there is also the 
question of time-scales and which non-linear couplings 
are more efficient than others. In experiments in David 
Weiss's group^ where there is currently no optical po- 
tential, a prethermalizcd GGE type state is found to per- 
sist for long times despite non-linearities arising due to 
effects such as band-curvature. This might be because 
the time-scales for the onset of dissipative and thermal 
effects due to these couplings are too long for experi- 
mental relevance. However we expect that application 
of a periodic potential will induce these effects more effi- 
ciently as a periodic potential is (dangerously) irrelevant. 
Once an experiment with an optical lattice is realized, a 
probe of the density-density response functions which di- 
rectly correspond to the correlators (e 10 ^" 101 '') that we 
have evaluated, should exhibit dissipative and thermal 
effects. 
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